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ABSTRACT 

The mean absolute extinction towards the central parsec of the Milky Way is 
Ak — 3 mag, including both foreground and Galactic center dust. Here we present 
a measurement of dust extinction within the Galactic old nuclear star cluster (NSC), 
based on combining differential extinctions of NSC stars with their vi proper motions 
along Galactic longitude. Extinction within the NSC preferentially affects stars at its 
far side, and because the NSC rotates, this causes higher extinctions for NSC stars with 
negative vi, as well as an asymmetry in the u/-histograms. We model these effects using 
an axisymmetric dynamical model of the NSC in combination with simple models for 
the dust distribution. Comparing the predicted asymmetry to data for ~ 7H00 stars 
in several NSC fields, we find that dust associated with the Galactic center mini-spiral 
with extinction Ak — 0.15 — 0.8 mag explains most of the data. The largest extinction 
Ak — 0.8 mag is found in the region of the Western arm of the mini-spiral. Comparing 
with total Ak determined from stellar colors, we determine the extinction in front of 
the NSC. Finally, we estimate that for a typical extinction of Ak — 0.4 the statistical 
parallax of the NSC changes by ^ 0.4%. 

Key words: galaxy center, nuclear cluster, kinematics and dynamics, dust, extinc¬ 
tion. 


1 INTRODUCTION 

Nuclear star clusters (NSC) are located at the centers of 
most spiral galaxies ( [Carollo et al.||1998| |Boker et al.|2002| ) . 
Their study became possible via high spatial resolution ob¬ 
servations from HST in the 1990s. They have properties sim¬ 
ilar to those of globular clusters although they are more com¬ 
pact, more massive and on average 4 mag brighter than the 


old globular clusters of the Milky Way (Boker et al. 2004 


Walcher et al. 20051. Many NSCs host an AGN (Seth et 


al. 20081 i.e. a supermassive black hole (SMBH) in their 


centers, have complex star formation histories (Rossa et 


al. 2006 Seth et al. 20061 and obey scaling-relations with 


host galaxy properties as do central SMBHs (Ferrarese et 
^[2^[Wehner fc Harris|2006[ |. 


The study of NSCs is of great interest because several 
of the most extreme physical phenomena occur within them 
such as SMBHs, active galactic nuclei, star-bursts and ex¬ 
treme stellar densities. The Galactic NSC is particularly in¬ 
teresting because of its proximity. At a distance of about 
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8kpc from Earth it is the only NSC in which individual stars 
can be resolved. 


The center of the galactic NSC harbors a SMBH (Gen- 


zel et al. 2010 Chez et al. 


Joint statistical analysis 


based on orbits around Sgr A* ( fCillessen et al.||2009| , star 
counts and kinematic data gives M, — (4.23±0.14) x 10 ®Mq 
and a statistical parallax Ro = 8.33±0.11 kpc ([Chatzopoulos 
et al.|20r5|). Recent studies (|Schodel et al.|2014||Chatzopou- 


los et al.|2015 1 have revealed that the NSC is flattened with 
an axial ratio q « 0.73, which is consistent with the kine¬ 
matic data ( [Chatzopoulos et al.||2015| . 

Observations of the NSC at optical-UV wavelengths are 
not feasible because of the high extinction Ay 30 mag due 


to interstellar dust (Scoville et al. 2003 Fritz et al. 20111. 


Therefore we rely on the infrared, with average K-band ex¬ 
tinction toward the central parsec close to Ak ~ 3 (Rieke 


et al. 1988 Schodel et al. 2010 Fritz et al. 20111. This is 


mostly foreground extinction; however, it is very difficult ob- 
servationally to measure the extinction variation along the 
line-of-sight within the NSC. 

The area around Sgr A* contains ionized gas which can 
be well described by a system of ionized streamers or fila¬ 
ments orbiting Sgr A* (Ekers et al.||198'3 Serabyn & Lacy 
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19851 that is also associated with hot dust |Lau et al.|2013 1. 


This complex structure of ionized gas is called the ’mini- 
spiral’ and consists of four main components: the northern 


arm, the eastern arm, the western arm and the bar (Zhao 
et al.|[2009l, surrounded by the circumnuclear disk of inner 


radius ~ 1.6pc ( Jackson et al.|19^ Christopher et al. 20051. 

In our previous dynamical study of the NSC, an asym¬ 
metry in the Wj-proper motions was observed in the his¬ 
tograms which was attributed to dust causing stars on the 
far side of the NSC to fall out of the sample. The aim of this 
paper is to present estimates for the extinction within the 
NSC based on our dynamical model and see to what extent 
this is correlated with the mini-spiral, to try to understand 
the slight asymmetry in the Vi velocity histograms of the 
NSC, and to check the impact of this on the fundamental 
parameters derived from the dynamical model such as the 
mass and the distance. 

In section 2 we discuss briefly our recent dynamical 
model of the NSC and describe qualitatively the effects of 
dust on the observed dynamics of the NSC. In section 3 
we show evidence based on the dynamics for the presence 
of dust within the NSC. In section 4 we develop a method 
for making an analytical model for the dust extinction that 
can be used on top of an existing dynamical model. Finally 
in section 5 we present extinction values for the dust within 
the NSC based on the prediction of the model in conjunction 
with the mini-spiral observations. 


2 EFFECTS OF DUST ON THE APPARENT 
DYNAMICS OF THE NSC 


In this section we give a brief description of the recent dy¬ 
namical model of the NSC from Chatzopoulos et al. (20151 
(Section 2.1) and we show initial evidence for dust extinction 
within the NSC (Section 2.2). 


2.1 Axisymmetric dynamical model of the NSC 


The dynamical model of the NSC described in |Chatzopoulos| 
et al. (20151 is an axisymmetric model including flattening 


and rotation which is an excellent match to histograms of 
proper motions and line-of-sight velocities in many bins on 
the sky. For the density model we used a two-component 


spheroidal 7 -model (Dehnen 1993 Tremaine et al. 19941 


which we fitted to the star count surface density data in {I, 


h) provided by Fritz et al. (20141. The inner rounder com¬ 


ponent can be considered as the NSC and the outer, more 
flattened component as the inner part of the nuclear stellar 
disk. Using the density we applied axisymmetric Jeans mod¬ 
eling in order to constrain the stellar mass M, , the black hole 
mass M,, and the distance Rq of the NSC which we found 
to be 

M.(r < 100") = (8.94±0.31| stat ±0.9| syst) X 10® M© 

M. = (3.86±0.14|,tat ±0.4|,y,t) x 10®M© 

Ro = 8.27±0.09| stat T 0.1|systkpc 

( 1 ) 

for the NSC only, not including the constraints from stellar 
orbits around Srg A*. 


Having this information we used the Qian et al. (19951 
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Figure 2. Figure illustrating the velocity profile asymmetries 
discussed in the text. Each small square represents a quadrant in 
the (Z, b) coordinate system, a ’+’ signifies an asymmetry in the 
vi VHs (e.g. right peak higher) a ’—’ the opposite asymmetry and 
’0’ no asymmetry at all. Only dust produces the same asymmetry 
in every quadrant. 


algorithm to calculate the even part of the 2 -Integral distri¬ 
bution function (DF) f{E, Lz). This allowed us to calculate 
the velocity profiles (VP) of the model. We found that the 
even part of the DF can predict very well the characteristic 
2-peak shape ( Schodel et al.|[2009 Fritz et al.||2014 1 of the 
velocity histograms (VH) for the Vi proper motion velocities, 
as well as the VHs in vi,. The addition of a suitable odd part 
in Lz to the even part of the DF represents the rotation of 
the cluster, so that also the uioa can be reproduced. 


2.2 Asymmetry of the vi proper motion 
histograms 

Upon a closer look at the velocity histograms in I direction 
(VH;) (see FigQ it is noticeable that the right peak is often 
slightly higher than the left peak. Seemingly there are more 
stars in the front of the cluster (positive velocities) than in 
the back. 

The proper motion data are from [lYitz et ah (20141. 
They are given in Galactic longitude I* and Galactic latitude 
b* angles centered on Sgr A*. In the following we always 
refer to the shifted coordinates but will omit the asterisks 
for simplicity. We assume that the rotation axis of the NSC is 
aligned with the rotation axis of the Milky Way disk. This 
is in accordance with the very symmetric Spitzer surface 
density distribution of [Schodel et al.| ( |2014[ ). 

At least three effects could produce an asymmetry as in 
Fig. Figure [^illustrates the effect of dust extinction with¬ 
out rotation, dust extinction with rotation, inclination and 
triaxiality on the VHs. Each small square represents a quad¬ 
rant of the shifted Galactic coordinate system {l,b) where 
Sgr A* is at the center. A ’-(-’ signifies an asymmetry (e.g. 
right peak higher), a ’ —’ the opposite asymmetry, and ’ 0 ’ 
no asymmetry at all (both peaks same height). We will see 
in the following Sections that dust with rotation produces 
the same asymmetry in every quadrant. Because of the dust 
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Figure 1. Histograms of u; proper motions and VPs predicted by an axisymmetric dynamical model of the NSC combined with a dust 
extinction model. The data and cells are as in |Chatzopoulos et al.|||2015[|. Each st ar is mapped to the first quadrant using {I, b) —>■ (|Z|, |b|). 
The smooth blue lines are derived from the recent model of |Chatzopoulos et al.| | |2015^ averaged over the cells combined as shown in the 
right panels, to which is added a homogeneous dust model with total Ak = 0.4, extending from —100'' to +100" along the line-of-sight(see 
Section 5). 


fewer stars will be visible at the back of the cluster, and 
because of the rotation the missing stars will be stars with 
negative velocities. Inclination of the NSC produces oppo¬ 
site results for the upper and lower quadrants, because when 
the line-of-sight does not pass exactly through the center, it 
passes through areas of unequal density in front of and be¬ 


hind the NSC. Finally, triaxiality produces opposite results 
for the right and left quadrants: consider a triaxial system 
with long axis in the Galactic plane but rotated away from 
the line-of-sight to the observer, with a triaxially symmetric 
rotation field. If the stars at the front of the cluster lead to 
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a higher peak for vi > 0 on the I > 0 side, say, then by 
symmetry the opposite will be true on the I < 0 side. 

Therefore, if we symmetrize the data to one quadrant 
only, the effect of the dust will remain while the effects of the 
inclination and triaxiality will cancel out, as illustrated in 
Fig-d This is what we observe for the symmetrized data of 
the NSC (Fig.Therefore we conclude that the observed 
asymmetry cannot be a result of inclination and triaxiality 
but that it may be produced by dust extinction in conjunc¬ 
tion with rotation. We note here that in order to observe 
this asymmetry (right peak higher than left) in the velocity 
profile in the I direction (VP;), the dust should be inside 
the cluster (i.e. within a few parsecs of the Galactic center) 
where the density is maximum, otherwise only a change in 
scale of the VPs would take place. 


3 DIFFERENTIAL EXTINCTION IN THE NSC 

For the work reported in the following, we use photometric 
data in the H and K bands and proper motions for 7101 
stars from Fritz et al. (20141. We split the data into a central 
and an extended field. The central field is a square centered 
on Sgr A* with size of 40” and contains 5847 stars. The rest 
of the stars belong to the extended fields, as shown in Fig. 

El 
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Figure 3. Map of Ax for the NSC derived from H — K colors. 
The map is similar to that of |Schodel et al.] | |2010[ |. The Galactic 
plane is shown as a black line. The central field is shown in white 
and the outer fields 01-03 in red, blue and brown. 


3.1 Total extinction 


We obtain the extinction towards each star from the H — K 
color (stars without H photometry and early-type stars are 
excluded). We obtain intrinsic color estimates by assuming 
that the stars are at the distance of the Galactic Center 
and that they are giants, as it is the case for most stars 


in the Galactic Center (Pfuhl et al. 20141. The intrinsic 


color varies between 0.065 and 0.34 but the majority of stars 
belong to a small magnitude range around the red clump. 
Therefore and also because the extinction is high, the in¬ 
fluence of intrinsic color uncertainties on the extinction is 
small compared to other effects, like photometric uncertain¬ 


ties. We use the extinction lawofFritz et al. (20111, implying 
Ak = 1.348F(i7-K). 

In this work we are mainly interested in the extinction 
variation Ak within the Galactic center. We obtain an es¬ 
timate for that by measuring for each star the extinction 
relative to its neighbors, more specifically relative to the 
median extinction of its 15 closest neighbors. Obvious fore¬ 


ground stars were already excluded in Fritz et al. (20141. 


By using 15 neighbors we obtain a robust median extinction 
estimate that is much less affected by extinction variations 
in the plane of the sky. To further reduce the influence of 
this extinction variation we exclude stars with too few close 
neighbors. 

Fig.i shows a map of interpolated total extinction for 
the central and extended fields based on H — K colors. The 
area within the white frame is the central field which is con¬ 


sistent with Fig. 6 of Schodel et al. (2010). Most of the ex¬ 


tinction of typically Ak = 3 mag shown in this plot is fore¬ 
ground extinction but a fraction ~ 0.4 mag or so is intrinsic 
to the NSC region as we show in the following. Fig.|^shows 
a histogram of the extinction for the central field. The mean 
extinction inferred from this plot is Ak = 2.94 mag with 
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Figure 4. Histogram of the extinction Ak based on H — K col¬ 
ors for all the stars. Mean extinction and standard deviation also 
given. Stars with small Ak are excluded because they are fore¬ 
ground stars ||Frit 2 et al.|2014ll. 


standard deviation 0.24 mag. This extinction value is very 
close to the value Ak = 2.74 measured by [Schodel et al.| 
(2010). The difference stems partly from the slightly differ¬ 
ent extinction laws and partly from the difference in areal 
coverage. Using the same law as Schodel et al. (20101 results 
in Ak = 2.85. 


3.2 Extinction in the NSC region 

The average differential extinction of stars as a function of 
Vi and Vb velocities is an important photometric quantity 
that can also be modeled and gives us information about 
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Figure 5. Average differential extinction of nuclear cluster stars 
plotted as a function of u;(blue) and Ui,(red) proper motion. The 
differential extinction is inferred from the difference in the color 
of a star to the median of its 15 nearest neighbors using the 
extinction law of [Fritz et al.| ( [2011^ and correcting also for the 
weak color variation with magnitude. 


the dust within the NSC. Fig. shows this for the central 
held. We notice that the average differential extinction as 
a function of vi is negative for positive velocities (preferen¬ 
tially at the front of the cluster) and positive for negative 
velocities (back of the cluster) i.e., stars at the front of the 
cluster are observed with less extinction than their neigh¬ 
bors. This hnding is consistent with the asymmetry of the 
VHs in I direction (Fig. and implies Ak — 0.4 within 
the NSC, see below. In contrast the average differential ex¬ 
tinction versus Vb is relatively hat and consistent with the 
symmetric bell shape of the VH in b direction. However we 
still notice a scatter of the points which is indicative of the 
systematic variations we should expect in Ak- 

Hot dust in the central parsec has been observed with 
mid-infrared emission (e.g., Lau et al. 20131. The emission 
traces dust at the inner edge of the circumnuclear disk 
(CND) and in the ionized streamers collectively called the 
minispiral (see Fig. 3, Lau et al.|2013 1. Fig. |^shows a vari¬ 
able extinction distribution in our central held which over¬ 
laps this region. This hgure was constructed by jScoville et al.j 
120031 using the ratio of their Paa emission measurements 
to the H92 q; radio recombination-line emission observed by 


Roberts & Goss (19931. The coloring signihes total extinc¬ 


tion including foreground along the line-of-sight and the con¬ 
tours show the outline of the mini-spiral. Dark blue colour 
indicates insufficient signal-to-noise in the H92a emission. 
The outline of the mini-spiral can be seen in this extinc¬ 
tion map, and the variability of the dust extinction along 
the minispiral suggests that a fraction of this extinction is 
likely to be associated with the mini-spiral itself and there¬ 
fore is located within the NSC. The outline of the minispiral 
in this map and in the MIR emission is closely similar, but 
the ratio of emission to extinction is not constant. In some 
regions outside the minispiral, particularly at larger radii, 
dust extinction is inferred from the ratio of Paa to 6cm ra¬ 
dio continuum emission (Scoville et al. 20031, presumably 
arising at least partially from colder dust in the surround¬ 
ing CND. These considerations will motivate our choice of 
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Figure 6. Extinction towards the ionized gas in the galactic cen¬ 
ter, as derived in jScoville et al.] ( |2003| their Fig. 5) from the ratio 
of their measured Pao emission to the H92a radio recombination¬ 
line emission of [Roberts fc Go^ ( [1993^ . The colours show ex¬ 
tinction in the V-band according to the scale shown; dark blue 
signifies no information because of insufficient H92a flux. The 
field is split into eight cells, taking into account the outline of the 
mini-spiral, the shape of the VHs in u; and the SAk curves. 


subfields for which we will construct separate dust models in 
Section For the development of the dust modeling in the 
next section we note that dust associated with the minispi¬ 
ral is likely to be concentrated in a small distance interval 
along the line-of-sight of order a fraction of the radius in 
the sky. 


4 DUST MODELING 

We saw in the previous section that the observed asymmetry 
of the VH in vi is likely to be associated with dust extinction. 
In this section we describe how one can make an analytical 
dust extinction model and use it with an already existing 
model of the NSC similar to that of section 2. 

For the rest of this work, along with I and b we use a 
Cartesian coordinate system {x, y, z) where z is parallel to 
the axis of rotation as before, y is along the line of sight 
(smaller values closer to the earth) and x is along the di¬ 
rection of negative longitude, with the center of the NSC 
located at the origin. 

First we need to model a luminosity function. We can 
do that by taking the product of two functions. The first rep¬ 
resents a power law function in luminosity, corresponding to 
an exponential magnitude distribution, L{m) = lO'*' "*. The 
second is an error function that represents the completeness 
function, so that: 


^ = L(m) X C{m) = 

107 >7i(i _ _ mo)/cr])/2 


( 2 ) 
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Here 7 is the power law index of the luminosity function and 
mo is the value where the completeness function C(m) has 
its half height. For the power law we set the index to 7 = 
0.27 ±0.02 as in Schodel et al. (2010). For the completeness 


function we set mo = 16.5 and a — 1 because we found that 
these values represent well the K luminosity data as shown 
in Figure The red curve of Figure shows equation 
with the chosen values. We have also investigated models in 
which we approximated the red clump bump at Kg ~ 16 in 
Fig. [7] by an additional Gaussian term. This led to ~ 10% 
increase in the derived Ak values. Since this is below the 
accuracy of our determinations, we use the simpler power- 
law luminosity function model in the following. 

Next we need the extinction variation over the line-of- 
sight which is just the derivative of the extinction over the 
line-of-sight i.e. daK/dy. The function daK/dy is general 
and could for example be represented as a sum of Gaussians 
but for simplicity we choose a square function, so that: 


daK 

dy 


0 , J/ < yi 
c,yi^y ^y2 
0 , 2 / > r /2 


(3) 


in which yi and 2/2 indicate the positions where the dust 
starts and ends respectively. The integral of eqj^over all line- 
of-sight is the maximum extinction Ak thus the constant c 
takes the value c = Ak /A 2 / where Ai/ = 2/2 — 2/i • Function 
[^integrates to: 


f 0 , y < 2/1 

o-K{y)^< ^ ( 2 /- 2 / 1 ), 2/1 s; 2 /^ 2/2 

I Ak , y > y2 


(4) 


p{y) = 


(5) 


The percentage reduction in observed stars as a function of 
line-of-sight distance is: 

f L{m — aK{y))C{m)dm 
f L{m)C{m)dm 

With this we can calculate the percentage reduction in num¬ 
bers of stars after extinction, Pmax = ^( 2 / 2 )- The percentage 
of stars hidden by extinction for Ak = 0.4 is about 25%. For 
the simple case where the luminosity function is a power law, 
the previous equation takes the form p{jj) = L{—aK{y))- 
The function for Ak = 0.4 is shown in figure 

Having that we calculate the VP; including the effect 
of dust extinction, denoted by VPD: 


VPD(i;i;a:,z) = - 


p{.y)ftot{E,L^)dviosdv^dy. ( 6 ) 


B>0 


Here E is the stellar surface density of the NSC model, and 
ftot{E, Lz) = fe{E, Lz) ± fo{E, Lz) is the total DF, consist¬ 
ing of an even part in Lz (contributing to the density) and 
an odd part (contributing to rotation), as in Chatzopoulos 
et al. ( 2015| ). Figurej^shows typical VPDs after adding dust 
with Ak = 0.4. We observe that the right peak of the VP; is 
now higher than the left peak which is a combined effect of 
dust and rotation. The dust does not produce an asymmetry 
for the VP in the h direction. One useful quantity is Ak(vi) 
averaged over the line-of-sight. This can be calculated with 


{Ak{vi-, X, z)) = 


Iff o-K{y)p{y)ftot {E,Lz)dVzdviosdy 

S X VPD (vi) 


(7) 


for each line-of-sight and connects our model with the pho¬ 
tometry. From this we calculate the average differential ex¬ 


tinction between stars of given vi and their neighbors, which 
corresponds to the data of Fig. 


6Ak{vi;x, z) = {Ak{vi)) - 


/(Ak (vi)}VPD (vi)dvi 
Ef VPD (vi )dvi 


( 8 ) 


In order to understand the effects of dust extinction on 
SAk{vi), we use two simple models for the dust distribution. 
The first is a homogeneous dust model that extends a few 
parsecs along the line-of-sight. The second is a thin ~ 10” 
screen of dust placed at several positions along the line-of- 
sight. We have verified that the width of the screen is not a 
sensitive parameter and the results are almost unchanged if 
we set it for example to ~ 5”. Figures]^ and 10 shows these 
effects. In Fig.j^we see the effect of three dust screen models 
and one homogeneous model on the VP; for Ak = 0.4. 

One important point to notice is that we can achieve 
the same effect (e.g. the same amount of asymmetry in the 
histograms) with several models by using different combina¬ 
tions of Ak and the distances at which the dust is placed 
along the line-of-sight, therefore the dust extinction model 
is degenerate. In Fig. we see plots of several screen and 
homogeneous models based on eq. The top panels show 
the SAk{vi) curves of a screen dust model placed at several 
distances in front of (left) and at the back (right) of the clus¬ 
ter for Ak ~ 0.4. The first thing to notice is that the further 
the screen of dust is placed from the center the smaller is 
the effect of dust. This makes sense since far from the center 
the density of stars is lower. We also note that if the screen 
of dust is placed in front of the cluster, the curves are close 
to constant for the stars behind the cluster (i.e., for nega¬ 
tive vi) since there is no dust there to affect the 5Ak{vi)- 
The opposite happens when the screen of dust is behind the 
cluster. The bottom left panel shows the shape of SAk{vi) 
for different Ak ■ The bottom right panel shows three homo¬ 
geneous models that extend over different distance intervals 
along the line-of-sight. In this case the curves are symmetric 
relative to zero. The dust extinction model was implemented 
with Wolfram Mathematica (|. 


5 PREDICTIONS OF THE DUST MODEL 

In the last section we described how one can include the 
effects of dust extinction in the dynamical modeling of the 
NSC, and calculated differential reddening signatures for the 
NSC stars. Here we proceed to model predictions and com¬ 
pare with both photometric and kinematic evidence. We will 
see that the asymmetries seen in the VPs for vi can be ex¬ 
plained as due to dust within the NSC. Our goal is to see 
whether dust in the Galactic Center mini-spiral can explain 
our data, and also to provide a rough extinction map of the 
central field based on the model and the available data. 


As discussed in Section 3.2 MIR observations (see the 
MIR map of jLau et al.|2013| ) show that hot dust is associated 
with the minispiral within the NSC although it cannot be 
assumed to trace the dust extinction directly. Figure [phased 
on Figure 5 of jScoville et al.| ( |2003| ) shows the total dust ex¬ 
tinction in the central 40” x 40” around Sgr A*, derived from 
the observed ratio of Paa to H92 q; radio recombination-line 
emission. This map includes both foreground and NSC ex¬ 
tinction. It shows the outline of the minispiral because it is 
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Figure 7. Modeled K-band luminosity function without and with 
dust extinction. The blue line shows a power law fit to the stars 
11 < K < 14 with power-law index of 0.27 di 0.02 as in |SchQdeI| 
|et al.| ( [2010| |. The red line shows the effect of the completeness 
function. The brown curve is the red curve shifted faintwards by 
= 0.4 dust extinction. The black line shows the K luminosity 
function for the Galactic center stars. 
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Figure 9. VPDs in I direction for different dust distribution mod¬ 
els, as indicated on the plot. The right peak is higher due to the 
combined effects of dust extinction and rotation. We can achieve 
the same effect (e.g. the same asymmetry) for several combina¬ 
tions of Aks and the distances where the dust is placed along the 
line-of-sight. All VPDs are for a line-of-sight with I = 10^^ and 
b = 10''. 
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Figure 8. Percentage reduction in observable stars after a dust 
effect with An = 0.4. as a function of line-of-sight for a homoge¬ 
neous dust model (red) and two dust screen models (brown, blue). 
About 75% of stars remain in the observed luminosity function 
after the extinction due to dust. 


limited by the signal-to-noise ratio (S/N) of the H92a recom¬ 
bination line flux, and the regions in which the extinction 
has an apparent value of Ay = 15 is where the S/N is not 
sufficient to derive the extinction. Fig. shows the greatest 
extinction in front of the Western arm, and relatively low 
extinction in front of the east-west bar. This suggests, in ad¬ 
dition to possible foreground extinction, position-dependent 
extinction associated with the minispiral itself. 

For our modelling, Fig.[^is split into eight cells and sub¬ 
cells. The reasoning behind the choice of these cells is based 
on the outline of the mini-spiral, but also on the statistics 
and shape of the velocity histograms in vi and the mean 
differential extinction variation 5Ak along line-of-sight in 
these cells. Some cells cover the region outside the minispiral 


where the observed ratio of Paa to 6 cm radio continuum 
emission ( Scoville et al.|2d03 1 also shows extinction, possibly 
in the foreground. This will also be tested by our modelling. 
The VHjs for all cells and the SAk as a function of vi and 
Vb are shown in Fig. for the following cells: 


• Cell Al: This cell’s small 5Ak values are consistent with 
the lack of features in the extinction map in comparison with 
other cells. Also the SAk values are consistent with the VH; 
since both peaks look symmetric which is a sign of lack (or 
small amount) of dust within the NSC. 

• Cell A2: The cell lacks strong mini-spiral features as cell 
Al. However both the VH; and SAk values show stronger 
effects of dust therefore this cell is separated from cell Al. 

• Cell B: The SAk values and the VHi are consistent 
with the strong features of the mini-spiral in the extinction 
map since the SAk points are higher and lower for negative 
and positive velocities respectively than the other cells and 
the asymmetry of the VH; is intense. We also note that 
the dust effects are similar within the whole area B because 
after splitting it into 2 sub-areas (not shown) we observed 
the same signature. 

• Cells Cl & C2: The cells Cl and C2 belong to the 
Northern Arm of the mini-spiral. The shape of the SAk 
data for both cells seems similar for negative velocities but 
the SAk for C2 is more symmetric and consistent with the 
extinction map hence we split the area into 2 halves. We 
notice also some asymmetry on the VH;s of both areas Cl 
and C2. 

• Cell D was separated from Cl & C2 because the SAk 
values look more symmetric than Cl & C2. We also note 
that the dust effects are similar within the whole area D 
because after splitting it into 2 sub-areas (not shown) we 
observed the same signature. 

• Cells E & F: The SAk values of these two cells look 
similar but the VH; of the cell F lacks the asymmetry char- 
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acteristic in contrast of cell E therefore we keep them sepa¬ 
rate. 


The central field with these eight cells is surrounded 
by a a more extended area with observations for about 2000 
stars. We split this area into three outer helds 01-03 placed 
around the central field as shown in Fig.[^ The VHj for these 
cells and the 5Ak as a function of vi are shown in Fig. 

Our goal is to give a model prediction of each of these 
cells (8-1-3 in total). For the model we use a thin screen of 
dust with width 10” because the dust associated with the 
minispiral is likely to be concentrated in a small distance in¬ 
terval along the line-of-sight. The precise width of the dust 
screen is not important as the dust signatures are insensitive 
to this parameter. The two main parameters are the total 
extinction in the screen and its location along the line-of- 
sight. However, as explained in the last section, these two 
parameters are partially degenerate. The degeneracy is par¬ 
ticularly strong for the VH; histograms. The 5Ak data in 
principle are sensitive to whether the extinction is in front 
or behind the Galactic center, as shown in Fig. |10| However, 
the Ak data have large scatter between adjacent data points 
such that points with seemingly small error bars can even 
have the ‘wrong’ sign of 5Ak (Fig. 111. This large scatter is 
also seen in the 5Ak versus Vb plots (also shown in Fig. 111 
where no dust signature is present. Therefore we decided to 
not try to fit the data using . 

Rather, we choose to place the dust screens along the 
line-of-sight according to other available information, and 
only deviate from this when this appears inconsistent with 
the shape of the 5Ak distribution. The total dust extinction 
of the dust screen is then chosen by eye mostly from the 
amplitude and shape of the 5Ak distribution, taking into 
account also the scatter of the 5Ak points, and to a lesser 
degree from the asymmetry of the VH; peaks. 

Specifically for the central field we use the three orbit- 
model of Zhao et al. (20091 for the three ionized gas struc¬ 
tures in the central 3pc (the Northern Arm, Eastern Arm, 
and Western Arc). We then map the center of each of the 
cells to a point in the relevant orbit plane according to its 
R.A and Dec. position. The distance from the center along 
the line-of-sight is given from the coordinates of that point 
on the orbit plane. For each cell, we use one common mean 
distance. Table shows to which orbital plane each cell is 
assigned, and the distance of the dust screen from Sgr A*. 

Fig.[TT] shows the predictions of the screen dust model 
for the 8 cells of the central held. The reasoning behind the 
choice of Ak and positions for the dust screen is based on 
the shape of the data and the examples of Fig. |10[ 


• Cell Al & B: We place the dust screen in front of the 
cluster according to the value of Table These two cells are 
interesting because they both show a correlation between 
the outline of the mini-spiral and the photometry data. They 
also exhibit the maximum contrast between the amount of 
dust. Cell B needs ~ 5 times more extinction than cell A. 

• Cell A2: The dust screen is placed behind the center 
according to Table 

• Cell Cl & D: The shape of the data in conjunction with 
the top right panel of Fig. [^indicate that the dust screen 
should be behind the cluster (also gives a much better 

in contrast with the value of Table A possible reason for 


Northern Arm 

E (-15”) 



Eastern Arm 

A2 (5”) 

Cl (-9”) 

C2 (-25”) D (-20”) 

Western Arc 

B (-25”) 

Al (-10”) 

F (-5”) 


Table 1. Each one of the eight cells of Fig. |10| belongs to an orbital 
plane | |Zhao et al.|200^ representing one of the three ionized gas 
(Northern Arm, Eastern Arm, and Western Arc) formations. The 
inferred line-of-sight distance of the dust screen from Sgr A* is 
given in the parentheses (negative points towards the earth). 


this could be that the single orbit description is not accurate 
for the centers of these cells. 

• Cell C2 & E: For both cells the dust screen is placed in 
front of the cluster according to Table 

• Cell F: The shape of the data in conjunction with the 
lower right panel of Fig. indicate that the dust screen 
should be centered. 


For the outer cells of the extended held the Ak is se¬ 
lected according to the general characteristics of plot |10| 


• Cell Of: the VHs and the photometry have a consistent 
signature. The extinction is close to Ak — 0.4 mag. 

• Cells 02 and 03: For these cells the photometry gives 
Ak — 0.3 and 0.35 but the peaks of the VHs are almost 
symmetric. 


Based on the model of the CND in Lau et al. (20131, we 
might have expected stronger extinction effects in helds Of 
and 03 than in 02. However, this is not conhrmed by Fig. |10| 
and the lack of correspondence between the 5Ak and the 
VHs is puzzling. 

Having the prediction for the cells we can estimate the 
foreground extinction for each cell using the Ak values of 
Fig.m The second row of table shows the total extinction 
of each cell of the central held based on Fig. The third 
row shows an estimate of the foreground extinction based 
on Ax (total) = Ax (foreground) -I- a; * Ax (NSC) where we 
approximate a; = 0.5 to take into account that of order half 
of the stars are unobscured because they are in front of the 
cluster. The exact value is not important for the scope of 
this paper but it should be closer to a; = 0.5 than a; = 1. 


6 DOES THE ADDITION OF DUST AFFECT 
THE MEASURED M., M* AND Ro? 

In this section we try to answer how much the dust within 


the NSC will affect the derived (Chatzopoulos et al. 2015 \ 


statistical parallax, supermassive black hole and stellar mass 
of the NSC. The foreground dust will affect the VPs only by 
a scale factor which does not impact the derived values. In 
[Chatzopoulos et al.| ( |2015| ) we derived new constraints on 
the Ro, M* and M. by fitting to the corresponding data the 
{v'^)\^bios the 2”“* order Jeans moments, that are 

moments of the even part of the corresponding VPs of the 
2-Integral distribution function. Therefore here the problem 
is reduced to how much the even part in Lz of the VPs 
changes after the addition of dust within the NSC. 

Figure |T^ shows the even parts of the VPs for vi,Vb 
and vios for the NSC dynamical model from jChatzopou-j 
los et al. (20151 with best mass and distance parameters. 
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Screen of dust Screen of dust 




i//(mas/y) 



-6 -4 -2 0 2 4 6 


t'/(mas/y) 



-6 -4 -2 0 2 4 6 


i;/(iiias/y) 


it/(mas/y) 


Figure 10. Shape of SAfc versus u; curve for several dust models. Top left: in front of the cluster for constant = 0.4. The numbers 
in the parenthesis show where the dust screen starts and ends respectively along the line-of-sight. Top right: Dust screen models placed 
behind the cluster for constant Afc = 0.4. Bottom left: Dust screen models placed slightly behind the cluster for several Ajy values. Bottom 
right: Homogeneous dust models that extend over different distances ranges along the line-of-sight for constant Ajy. The line-of-sight 
for all curves has coordinates I = 10" and b = 10". 


Cell 

A1 

A2 

B 

Cl 

C2 

D 

E 

F 

Total extinction 

2.73 

2.99 

3.1 

2.82 

2.87 

2.78 

2.84 

2.95 

Foreground extinction 

2.66 

2.79 

2.7 

2.67 

2.72 

2.61 

2.69 

2.8 


Table 2. Extinction values per cell based on Fig. The foreground extinction of each cell is estimated according to Ajy (total) = 
Aif (foreground) -ha:* Ajf (NSC) where we approximate x = 0.5 to take into account that of order half of the stars are unobscured because 
in front of the cluster. The exact value is not important for the scope of this paper but it should be closer to a; = 0.5 than a; = 1. 


and for the same dynamical model including the screen dust 
prediction, for a line-of-sight through cell B which has the 
largest amount of extinction {Ak = 0.8) among the cells of 
the central and outer fields. We notice that the difference 
between the VPs for this amount of extinction is very small. 
Specifically the average difference of the 2"'^ moments of the 
VPs between the two models is ~ 1.5%. If instead we use 
Ak = 0.4, close to the average extinction within the NSC 
inferred from this work, the difference is smaller than 0.5%. 
The relative differences of aiaajo^, and ctjos/ct; between the 
model with no dust and the model with Ak = 0.4 are simi¬ 


larly small, 0.2% and 0.6%, respectively. Therefore we con¬ 
clude that the systematic effects on the statistical parallax 
due to dust are within the estimated errors of|Chatzopoulos| 
[eFUIpoTsl ), causing the distance to the NSC to decrease 
by ~ 0.4% ~ 30pc. 

That the changes in the even part are so small can 
be explained by the following formal argument for the 
VP;. We show that for small amounts (1®* assumption) of 
homogenized (2"''^ assumption) dust around the center the 
even part of the renormalized VP in vi is the same as that 
for no dust. The odd part is a direct indicator of dust at 
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i'/(iiMsy) i'i(mas/y) 


Figure 11. Predictions of the model with the VH and 5Ai^ data for each cell. The numbers in the brackets show where the screen of 
dust is placed relative to the center. The line-of-sight is placed at the center of each cell. Reduced a^lso provided for the histograms 

and the photometry. For the cells Al, A2, B, C2, E, the screen dust distance is based on the orbit models of the mini-spiral (|Zhao et al.| 
[ 200 ^ . 
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I'l (iiias/y) t';(iiia,s/y) 




1'; (iiias/y) 


t-,(nias/y) 


Figure 12. Predictions of the model with the VH and SApc data for the outer cells shown in FigUl The numbers in the brackets show 
where the screen of dust is placed relative to the center. The line-of-sight is placed at the center of each cell. Reduced 
provided. 


that point (/, b) and can be used to estimate Ak^ 

When the luminosity function is a power law and the 
dust is homogeneously placed i.e. yi = —y 2 in equations 
and[^ p{y) takes the form within the dust area: 

I ^Ky \ 

p{y) = L{-aK{y)) = 10 2 Ay j ( 9 ^ 

For small Ak <S ( 7 ln( 10 )) ^ ~ 1.61 we have: 


Therefor^the function g{y) = p{y) — p(0) is odd every¬ 
where (including the area where there is no dust). Next for 
simplicity we use ftot{E,L^) = f{E,xvios - yvi) -A f{viy) 
because y appears within / only with the form of y^ and 
viy. 


p{y) ~ 1 - 7 In(lO) ^ - 7 In(lO) ^j/ 

= p( 0 ) - ky 

where p(0) = 1 — 7 ln( 10 )^^ and k = 7 ln( 10 )^^. 


^ we showed this for the case where the luminosity function is 
a power law but the same holds in the general case Lira) where 

7 ■ In(lO) is replaced by ■ 
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Figure 13. Even part of the renormalized VPs for vi,vi) and vios 
for the NSC dynamical model with no dust from |Chatzopoulos| 
|et al.| | [20T^ , and for the same dynamical model including the 
dust screen prediction based on cell B. The black line in the up¬ 
per panel corresponds to the odd part of the VP^ The VPs are 
calculated at the center of cell B. 


Now we have: 

1 1Piy) II ifiviy) + f{-viy)) dv^dviosdy 
-IIp(O) II (f(viy) + f{-viy)) dv^dviosdy = ( 11 ) 

1 19{y) II ifi^iy) + fi-viy)) dv^dviosdy = 0 

The previous is 0 because g{y) is an odd function and 
II{fi'^iy) + f{~'^iy))dVzdvios is an even function of y 
therefore: 

1 1 Piy) II ifiviy) + fi-viy)) dv^dviosdy = 

( 12 ) 

ip(0) III ifiviy) + fi-viy)) dv^dviosdy 


And thus: 


VPDeven(U!) =p(0)VPl(Ui) 


(13) 


To hnd the constant p(0) we integrate once more over 
the velocity this time: 

j VPDeven(U!)dUi = p(0) J YPl{vi)dVl = p(0) (14) 

Therefore p(0) is the normalization factor of the VPDi. The 
above means that to first order, dust does not affect the even 
part of the VPD in I direction significantly except for a scale 
factor. Fig[T3] shows that the effects for VPDb and VPDios 
are similarly small. If extinction within a stellar system is 
small enough {Ak <C 1.6 for the NSC) then fitting the even 
part of a model’s VPs to the even part of the VHs is sufficient 
to get accurate estimates of the M., M* and Ro parameters. 


In Chatzopoulos et al. (2015 I we used the root mean square 


velocities that are moments of the even parts of the VHs 
to 6t the M,, M* and Ro parameters of the axisymmetric 
model. Therefore we expect that their values will not be 
affected by dust more than 0.4% as we explained earlier. 

In principle the odd part of the VPDi can be fitted to 
the odd part of the VHs and this part is scale free since the 
scaling factor is already known from the even part therefore 
one can fit the Ak for some combination of cells. 


7 DISCUSSION & CONCLUSIONS 


The main goal of this work was to understand the slight 
asymmetries in the VPiS of the NSC and their influence on 
the dynamical modeling following the recent work of [Chat-] 
[zopoulos et al.| ( |2015| ). Our interest was triggered by the 
observation that the right peak of the VH; is often slightly 
higher than the left. A plausible explanation was given based 
on the existence of dust within the NSC. Because of the dust, 
the apparent number of stars behind the NSC is smaller than 
that in front of the cluster. This in conjunction with the ro¬ 
tation can explain the observed characteristic. 

In order to quantify the dust effects, we worked with 
proper motions and photometry for ~ 7100 stars from |Fri^ 


et al. (20141. We applied an analytic dust extinction model 


together with our recent NSC dynamical model. The extinc¬ 
tion model gave us reasonable results and was able to predict 
both the signature in the VPs and the photometry. 

Observation of the NSC in the optical is almost impos¬ 
sible because of 30 mag extinction. In the infrared the situ¬ 
ation is much better since Ak ~ 3 mag. Most of this extinc¬ 
tion belongs to the foreground and does not affect the shape 
of the VHs (except a normalization factor). We find here 
that a small fraction of the total extinction value (~ 15%) 
belongs within the NSC. 

The area between ~ 1 — 1.5pc radius consists of several 
streamers of dust, ionized and atomic gas with temperatures 
between 100A — 10"^A and is called ’’ionized central cavity” 

( Ekers et al.jl983 1. The mini-spiral is a feature of the ionized 
cavity, and is formed from several streamers of gas and dust 
infalling from the inner part of the CND ([Kunneriath et ah] 


20121. It consists of four main components: the Northern 


arm, the Eastern arm, the Western arm and the Bar (Zhao 
|et al.|2009| that can be described well by streams of ionized 


gas or filaments orbiting Sgr A* (Serabyn fc Lacy|1985l. 
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We investigated how well the mini-spiral correlates with 
the extinction effects in the NSC data within the central 
field. To assess this we first investigated whether our extinc¬ 
tion model puts the dust on the same side of the NSC as does 
the mini-spiral interpretation. This is true for six out of the 
eight cells (except cell C2). We found the largest extinction 
{Ak = 0.8) in cell B where also the largest extinction is in¬ 


ferred from the extinction map of Scoville et al. (20031, and 


particularly low extinction in cell A where the extinction 
map is consistent with only foreground extinction. 

We can estimate the mass that corresponds to a given 
amount of extinction using 


Pd 


Ak 


Ay 1.086ka 


(15) 


where pd is the dust density, ka = 1670cm^/g (Draine 
2003a I is the mass extinction coefficient for the K-Band, 
Ak = 0.4 from the model prediction and Ay = 10” is the 
width of the dust screen. We find that ~ 1.8x 10“^®kg/m® 
and the dust mass within a parallelepiped with dimensions 
(10”, 40”, 40”) centered on Sgr A* is Md ~ 3Mq. Since the 
dust extinction model presented here is not precise we con¬ 
sider that this estimate is correct only within an order of 
magnitude. This value is within the range of 0.25 — AMq for 
the mini-spiral found from other works | |Zylka et al.||1995{ 
Kunneriath et al.|2012 Etxaluze|2011 1. 

Finally we showed that for small values of extinction 
the even parts of the VPs are not affected significantly. 
As a result, the measured M, and M* parameters of 


Chatzopoulos et al. (20151 do not change by more than 
~ 0.4% for extinction Ak — 0.4 and the inferred Rq is 
decreased by the same amount, which is less than the 
smallest systematic error (for the statistical parallax) 
inferred for these parameters in Chatzopoulos et al. (20151. 


Our results can be summarized as follows: 

• We showed that extinction due to dust explains kine¬ 
matic asymmetries and differential photometry of the NSC, 
and measured the amount of extinction within the NSC by 
combining a dynamical model with a dust extinction model. 

• We presented an extinction table for the dust within 
the NSC in several cells. 

• We found that the distribution of the dust is consistent 
with the extinction being associated with the mini-spiral for 
six out of eight cells. 

• Systematic effects due to dust with typical extinction 
Ak — 0.4 affect the M., M* and parameters deduced 
from previous dynamical modeling only by ~ 0.4%, which 
is smaller than their estimated systematic errors. 
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